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Abstract 



We study the ensemble-averaged conductance as a function of applied magnetic field for ballistic electron 
transport across few-channel microstructures constructed in the shape of classically chaotic billiards. We 
analyse the results of recent experiments, which show suppression of weak localization due to magnetic 
field, in the framework of random-matrix theory. By analysing a random-matrix Hamiltonian for the 
billiard-lead system with the aid of Landauer's formula and Efetov's supersymmetry technique, we derive 
a universal expression for the weak-localization contribution to the mean conductance that depends only 
on the number of channels and the magnetic flux. We consequently gain a theoretical understanding of 
the continuous crossover from orthogonal symmetry to unitary symmetry arising from the violation of 
time-reversal invariance for generic chaotic systems. 
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1 Introduction 



Random- matrix theory has proven to be a very successful tool for the understanding of localization phe- 
nomena in disordered diffusive systems, for example, universal conductance fluctuations H, Altshuler- 
Aronov-Spivak oscillations ^, Q and persistent currents In a number of recent experiments, weak- 
localization effects have been observed in ballistic cavities in which electron scattering from the 
boundaries is classically chaotic. It is interesting to also investigate the application of random-matrix 
theory to such problems as one can expect it to be able to address questions not accessible to other 
methods, such as semiclassical analysis. 

In this work, we look at ballistic electron transport through a microstructure that has the shape 
of a classically chaotic billiard, in the presence of an applied magnetic field B. The typical set-up for a 
stadium is shown in Figure 1. The stadium, whose length L is less than the mean free path £ and the 
dephasing length L^, has attached to it two asymmetrically positioned ideal quasi-one-dimensional leads 
which carry the incoming and outgoing electron transverse modes. We assume an equal number, M/2, 
of incoming channels (labelled by index a) and outgoing channels (labelled by index b). Thus M counts 
the total number of open channels coupled to the stadium. 

We use this physical system to study the breaking of time-reversal symmetry, as measured by the 
conductance G{B) = {e^ / h)g{B), as the magnetic field B is turned on. The relevant experimental ob- 
servation |6|-| is that the average dimensionless conductance g dips at zero magnetic field B = 0. The 
g versus B curve is illustrated schematically in Figure 2. The conductance interpolates between a min- 
imum value, which can be associated with the Gaussian orthogonal ensemble (GOE), at _B = and an 
asymptotic upper value, which can be associated with the Gaussian unitary ensemble (GUE), for large 
B. The difference Sg between the upper GUE limit and the mean dimensionless conductance g represents 
the weak-localization contribution to g. The suppression of weak localization due to non-zero magnetic 
field, i.e. the fact that Sg a.s B oo, corresponds to a GOE GUE transition. 

The appearance of the orthogonal and unitary random-matrix ensembles in this discussion arises from 
the fact that the fluctuation properties of the quantum spectra for closed chaotic billiards are the same as 
those of (i) the GOE for B = 0, and (ii) the GUE for B sufficiently different from zero. The former, being 
an ensemble of Gaussion distributed real-symmetric matrices, respects time-reversal invariance, while the 
latter, being an ensemble of Gaussian distributed Hcrmitian matrices, does not. 

The presence and magnitude of the weak-localization term in 'g{B) at i? = can easily be understood 
in the framework of random-matrix theory (without detailed calculation) as a consequence of (i) elastic 
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enhancement (or equivalently, coherent back-scattering), derived in Section 5.1, and (ii) imitarity of the 
S-matrix. Our starting point is the connection between the conductance g and the S-matrix given by 
Landauer's formula: 

M/2 M 

E {\Sat\' + \Sta\'} , (1.1) 

a=l b=M/2+l 

where a factor of two for spin is effectively included in this expression. Here, the index a sums over the 



M/2 incoming channels and b sums over the M/2 outgoing channels, so that the summation in Eq. (1.1) 



runs over two non-overlapping sets of channels {a} =/= {b}. If we assume all channels to be equivalent, so 



that all mean values jS'abp are equal for a ^ b, then we obtain 



M\ f M\ ——^ M^ 



9-'^\'-f) [f)\SabV^—\Sa,V- (1-2) 



Now, according to (i), we have 



with 



As a consequence of (ii), we have 



\S,cV = n\Scc'V for c^c' (1.3) 



2 for GOE 
1 for GUE. 



(1.4) 



M 



J2 l^cc'P = trlM =M. (1.5) 



c,c' — 1 

Combination of the two equations ( |1.3| ) and ( |1.5| ) yields 

M M 

C=l c,c' = l 



MK\Sab\^+M{M~l)\Sab\'' (1.6) 



for any a ^ b. Therefore, we see that 



1 _ 1 1/M for GUE (^^ = 1) 

~ M -l + n^\ l/{M + l) for GOE ^ 2) . ^^-'^ 



We can now use Eq. (1.7) to look at the expansions of the corresponding values of g in powers of Af ^. 



For the pure GUE (large magnetic field), we obtain 

5GUE = M/2 , (1.8) 
while for the pure GOE [B — Q), we see that 

5goe 



2(M + 1) 

= ^cue4 + o(^) . (1.9) 

The contribution —1/2 to the RHS of the equation above represents the weak- localization term in the 
conductance, and gives the magnitude of the gap between the upper and lower asymptotic values in 



Figure 2 in the limit of large channel number M . Weak localization is absent for the GUE, evidently 
being suppressed by the magnetic field. 

One should also observe that, for general values of Ad , the weak localization term has the value 
M/{2{M + 1)). It increases from 1/3 for M = 2 to 1/2 for M — > oo. Nonetheless, the weak localization 
correction decreases relative to g which is linear in M. All this agrees with recent results obtained |^ for 
the circular unitary (CUE) and circular orthogonal (COE) ensembles. The agreement is not surprising: 
It has been shown |l^ that the GOE and the COE distributions for the S'-matrix elements coincide when 
the coupling with the leads is maximized. Similar arguments apply to the unitary case. 

We now consider the question: How does a finite magnetic field interpolate between 5qoe (-^ ^ 0) ^^^^ 
5gue (-^ sufficiently large)? Semiclassical analysis [|l^ implies that Sg{B) — ^gue ~ Vi^) ^ Lorentzian. 
However, this work was only able to include the diagonal contributions (i.e. taking into account only in- 
terference between symmetry-related classical paths) , and was not able to achieve quantitative agreement 
with fully quantal numerical simulations. In view of this, one should ask whether the full result remains 
Lorentzian. 

Our aim in this investigation is to understand the generic features of the continuous GOE — )■ GUE 
crossover as a function of magnetic field B for any classically chaotic microstructure in the framework 
of random-matrix theory. From the semiclassical perspective, our random-matrix approach may be 
viewed as taking into account the generic properties described by the long classical trajectories, in which 
system-specific information is washed out. Therefore, we shall need to restrict attention to billiard/lead 
configurations that are dominated by long classical trajectories. Our results will not take account of 
effects associated with the specific geometry of a billiard, which would be probed by short trajectories. 
Moreover, the density of long trajectories is large — growing exponentially with the period — and 
therefore virtually impossible to account for in the framework of semiclassical theory. 

On an energy (rather than a time) scale, such generic features are displayed on intervals measured 
in units of the mean level spacing d. The suppression of weak localization becomes effective whenever 
the matrix elements of the interaction of an electron with the magnetic field become comparable with 
d. This is why we believe that our approach should yield results which are relevant for the problem. 
There exists a time scale — the mixing time T,„ix, defined below — which relates to the strength of these 
matrix elements. The suppression of weak localization will be seen to depend actually on the competition 
between this time scale and a second time scale Tdoc, the mean life-time of the levels in the billiard. The 
associated width Fdcc — ft/Tdcc grows linearly with the number of channels. Our final result essentially 
depends on the ratio of these two time scales. 
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The paper is set out as follows: In Section 2, we define the model and discuss the values of its 
parameters. Section 3 deals with relating the parameter that controls the GOE GUE crossover to the 
strength of the external magnetic field. In Section 4, we outline the supersymmetry formalism used to 
extract an exact result for the ensemble-averaged conductance, and explicitly construct the generating 
function that is central to this approach. We also summarize the steps leading to the final analytical 
expression for g, which is then displayed and discussed in Section 5. Results of the numerical evaluation 
of the magnetic field and channel number dependence of the weak-localization term are presented in 
Section 6, along with a discussion of the Lorentzian approximation. Section 7 contains some interpretation 
of our results, as well as a comparison with experimental data. Finally, we present our conclusions in 
Section 8. A few appendices contain various technical details. 



2 The Model 

We first consider an electron moving in a closed billiard, whose Hamiltonian H we model in terms of an 
ensemble of random matrices. We then take account of the presence of external channels by introducing 
a coupling between H and the leads. This allows us to define a scattering matrix which we relate to the 



conductance via Landauer's formula, viz. Eq. (1.1). We shall maximize the coupling between H and the 
leads. Then, the (ensemble-averaged) conductance g(t, M) depends only on the strength of the magnetic 
field B (expressed in terms of a dimensionless parameter t), and on the total number M of channels. 
Using ergodicity, we compare this quantity with experiment. We neglect the Coulomb interaction between 
electrons as well as all other dephasing effects. 

For the Hamiltonian of the closed billiard, we consider an ensemble of random matrices H^i, of 



dimension N, given by p4| 



H,.=Hl^J+^J-LHlt^. (2.1) 



The ensemble _ff^^^ represents a GOE with zero mean and whose non-diagonal elements have variance 
X^/N, while i?^^^ denotes an ensemble of Gaussian distributed real antisymmetric matrices whose non- 
vanishing elements have zero mean, the same variance /N and are uncorrelated for fi > v. Explicitly, 



H^f^^H^^l, - —{5^_,^'5yy' - 5^y>5^^') . (2.2) 

The parameter A fixes the mean level spacing d of the GOE and must, in principle, be adjusted to the 
local mean level spacing at the Fermi energy of the billiard at hand. In the end, we shall take the limit 
N CO, in which case the average level density of the ensemble (^.l[) will take the form of Wigner's 
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semicircle. We shall also normalize all energies such that the Fermi energy is zero, Ep = 0, in which case 
it sits at the centre of the Wigner semicircle. The mean level spacing at the Fermi energy is then given 
by d = nX/N. 



To estimate the scale of the GOE — ^ GUE crossover, we note that this transition will take place 
when the typical matrix element of the time-reversal violating perturbation becomes comparable with 
the mean GOE level spacing: 

which implies t ^ 1, i.e. for a value of t of order unity on the scale of N. This is why we chose to include 



a factor N-^^'^ in front of H^^^ in Eq. (O) 



The coefficient ofiH^^^ gives the strength of the time-reversal violating perturbation and is expected 
to be proportional to the magnetic field, since the variation of the levels of the stadium with magnetic field 
is linear in the flux cj) = BA through the area A of the stadium. It is convenient to express \/t = k(j)/(j)o, 
where (po — h/eia the fundamental flux quantum, and k is then a number of order unity. An estimate 
of k can be provided by computing the As-statistic of a closed stadium for various values of the applied 



magnetic field, and then fitting the results using the matrix ensemble H^i, of Eq. (2.1). More details of 
this procedure will be given later. 

Next, we must consider the coupling of the closed stadium with two ideal external leads. The 
Hamiltonian operator for the total system (comprising stadium and leads) is given by 

N M 



+ ME, c) {\E, c)W,^{E){^l\ + h.c.} . (2.4) 

The first term represents the stadium Hamiltonian as discussed above. The second term describes the M 
free electron channels in both of the leads. Each channel c is defined by a propagating transverse mode 
with an energy E below the Fermi energy. Thus \E,c) is the plane wave mode of transverse energy £c 
and longitudinal momentum hk, such that the total energy of the state is E — ec + TT?k'^/{2m*), with m* 
being the effective mass of the electrons, and dfi{E, c) is the corresponding energy measure: 



d^i{E,c) = J-!I^Q{E-e,)^^= , (2.5) 

V STT^n \J E — €c 

where Q{E) denotes the step function. The final term contains the lead-stadium couplings. We have 
denoted by Wc^i the matrix element which connects the basis state in the stadium's internal region 
with channel c, where /i = 1, 2, . . . , and c— 1,2,..., M. We assume that the breaking of time-reversal 
symmetry by the magnetic field B need only be taken into account in the Hamiltonian (^.l|) of the stadium 



and not in the matrix elements Wc^ describing the passage of the electron between one of the channels 
and the stadium interior. Accordingly, we set Wcf^ = W*^. 

General quantum scattering theory allows us to write down the scattering matrix Scd for the problem 
at hand It has the form 

Scd - Scd - 2tTrJ2Wc,.[D-\,Wd. , (2.6) 
where the inverse propagator D^^ is given by 

J2 Wc^cWc, . (2.7) 

c 

Here, the energy E is taken to be the Fermi energy Ep, and thus will ultimately be set to zero. The 
scattering matrix defined in Eqs. (p^GpJ) is manifestly unitary but not symmetric (unless t — 0). In 



writing eq. (2.7), we have assumed that the energy dependence of the matrix elements Wcf^ is negligible 
over an energy interval of the order of d, and we have therefore omitted a principal-value integral. It 
was shown in Ref. |l^ that taking such an integral along complicates the calculation but does not affect 
the GOE scattering result, except for a redefinition of the transmission coefficients introduced below. 



Insertion of Eqs. ( p.(]| ) and (2.7) into the Landauer formula (1.1) furnishes us with a formal expression 
for the conductance g in terms of the Hamiltonian, which must now undergo ensemble averaging. For 
this purpose, it is useful to note the compact representation 

SabS:, = ATrn'^D-'n'^iD-']^ , (2.8) 

where the trace 'Tr' runs over the level index /i and we have introduced the quantities 

f^M- = T^Wc^Wc, . (2.9) 

It is worth pointing out at this stage that the parameters {Wc/^} will not appear explicitly in our final 
expression for g, for the following reason: The random-matrix ensemble ( |2.lD is invariant under unitary 
transformations. Therefore, only unitary invariants constructed from the Wcf^^s can appear in averages 
like The only such invariants are the bilinear forms 

^cd^^^'E^^M^'^/^- (2.10) 

Because of the summations in Eq. (^]|), we may assume without loss of generality that 

Vcd = (2.11) 

is diagonal in the channel indices if c and d are channels in the same lead. In assuming that this relation 
holds in general, we disregard short classical trajectories connecting the two leads. It turns out ||lj] that 



the VcS enter into the final expression for g only in the form of the dimensionless coefficients 



(2.12) 



(A + X,)2' 

where Xc — irNv'^. Likewise, the parameter A specifying the local level density appears only in the 
combination Xc- Obviously, we have < Tc < I for all c. The quantities Tc are interpreted as transmission 
coefficients which measure the strength of the coupling between the cavity and the leads. The absence of 
barriers between cavity and leads motivates us to put Tc — 1 for all c, thereby maximizing the coupling 
to the leads. Then, g depends only on B and on M . All dependence on A and on the Fermi energy has 
disappeared, and we are left with a calculation that contains no free parameters. 

We emphasize that in our approach, all channels are strictly equivalent. There is no distinction made 
between channels in lead one and those in lead two. Hence, the sums in Eq. (^^) are restricted only in 
the sense that they must be carried out over different sets {a} ^ {b} of channels. Clearly, this cquivalenve 
assumption is justified only for long delay times of the electron inside the billiard. The experiments [^|| 
are arranged in such a way that it is unlikely for electrons to pass directly from one lead to the other, 
lending credence to this assumption. 

A crucial constraint applies in our calculation of 5. We cannot use the usual perturbative expansion 
in powers of 1/M as done in previous work First, experiments typically have a small number of 
open channels, Af = 2,4 or 6. For example, in the experiments of Refs. ^ and 0, the stadium billiard is 
connected to two almost one-dimensional contacts in each of which the number of channels M /2 depends 
on the gate voltage and varies between one and three. For Af 3> 1, the semiclassical approximation is 
expected to work. For small M , on the other hand, we enter a different regime which is not accessible 
to perturbative methods. Physically speaking, the levels in the billiard are broadened by the coupling to 
the leads. For small M , this broadening is of the order of the mean level spacing d. This is the regime we 
mainly address in this paper, although our results are generic and apply to all values of M. Second, and 
more importantly, we shall see that the dependence of the mean conductance enters in the combination 
t/M. This implies that the perturbative results will be good only as long a,s t M, but this will be 
shown to fall short of the crossover region. Consequently, we must resort to a computation that is exact in 
the number of channels M . This we shall achieve by using the method of the supersymmetric generating 
function [|l5[ Q to find an expression for g(t, M) valid for all values of t and M. 
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3 Magnetic- Field Parameter 



In order to completely specify our model in terms of physical quantities, it remains to discuss how one can 
relate t to the strength of the magnetic field B. A previous investigation of the GOE GUE transiton 
in terms of the same magnetic-field parameter t by Pandey & Mehta jl6| offers the possibilty of relating t 
and B by direct inspection of the spectral fluctuations of the stadium for different values of the external 
magnetic field. Also, in more recent work, Ozorio de Almeida |l7|| derived semiclassically the dependence 
of t on i? for a chaotic system in terms of purely classical correlators. Although these classical quantities 
are difficult to compute with very good accuracy, an efficient average over the phase space is sufficient to 
provide another independent quantitative estimate of t in terms of physical quantities. In this section, 
both approaches are discussed. 

We shall start with the first method, and consider the Schrodinger equation for an electron mov- 
ing in a billiard in the presence of an external magnetic field B. Working in the gauge A{x,y) = 
{—By/2,Bx/2,0) and introducing rescaled coordinates x' — rj^^^x, y' = ij^/^y with t] = in/ A, A being 
the area of the stadium, we arrive at the Schrodinger equation in the form 



/ \ 2 

2 i (j) ( I d I 9 \ TT ( <j> \ , ,2 ,2s A 2m 



^P{x\y')^0, (3.1) 



where — h/e is the magnetic flux quantum. The boundary conditions are given hy ip{x' ,y') ^ for 
{x',y') £ dD, where dD is the border of the stadium. This is the eigenvalue equation that we solve 
numerically. We subsequently drop the primes on the coordinates. Recalling the Weyl formula for the 
average level density p^^ for Dirichlet boundary conditions in the absence of magnetic field. 



A 2m r 2m 



where V denotes the perimeter of the stadium, we see that, up to a correction of order E the rescaling 
leads to to an energy spectrum of average density one in units of 2m/ h. 



The linear term in (j) in Eq. (3.1) breaks time-reversal symmetry. The term quadratic in brings 



the system to the integrable limit as one increases cf). This term is of little importance for the range of 



B fields relevant to the experiment under consideration. The structure of Eq. (3.1) already suggests a 
linear dependence between the parameter t and (t>/(j)o, provided that the operator O = i(xd/ dy — yd/dx) 
locally connects different states efficiently, thereby mimicking a random matrix. Although the latter 
assumption seems to be rather natural, in a dynamical system, the variance of the matrix elements 
associated with O will typically vary with the energy. A qualitative prediction of the trend is possible by 
invoking the following arguments: For B = 0, Berry conjectured that the wavefunction of an eigenstate 
of energy E — fi^k'^ /2m for a classically chaotic billiard is typically given by a superposition of plane 
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waves a.i exp(iki-r) with randomly oriented, k= jk^j and Gaussian amplitudes Oj. Since according 
to the previous conjecture, O — —h{r x k)^, it follows that for eigenstates and v with k,y ^ k^^ ^ k, 
O ~ — ?ifc(i/|(r X k)z|/i) = kv^fi, where k = k/fc and is a random variable. In order to compute the 
wavefunctions of the stadium with good accuracy. Heller typically uses kV plane-wave components. 
Taking into consideration the number of components and the randomness of a^, one can estimate that 
Vp,i, cx k^^. From this naive argument, one concludes that oc k. Therefore, any numerical study of 
the GOE GUE transition in a dynamical system will be quantitative within a given energy range. A 
precise correction for the energy dependence of our results can be made via the semiclassical approach, 
and will be discussed later. 

Analytical results for the GOE GUE transition are available for the fluctuating part of the level- 
level density correlator, defined as 



Y2{e/d) = d^p^{E)p^iE + e) (3.3) 



where d is the mean level density. For the model defined in Eq. (2J^), one finds in Ref. |l^ the following 
analytical expression for Y2 as a function of x = e/d and t: 

Y,i,,t) ^ - \ fdrre-^V^^in™ T ^e-^'^^^ . (3.4) 

Unfortunately, this statistical measure is not very efficient for quantifying t in practice. Many numerical 
studies show that a more adequate one is the A3 (a;) statistic, which can be related to l2(x) straight- 
forwardly We have proceeded as follows: For different values of 0/00, we computed the A3 level 
statistics and adjusted t in Eq. ( |3.4[ ) to obtain the best fit. The results for the quarter stadium, with 
radius R = 1 and length £ = 0.5 are presented in Figs. 3a-3c. The statistical analysis of the spectra 
was performed over an energy interval where the level density stays almost constant. (Only the leading 



term in Eq. (3.2) is relevant.) Since we are interested in a quantitative determination of the transition 
parameter, care must be taken with respect to the following: (a) Berry's saturation of the A-statistics, 
which would jeopardize the method. By restricting our analysis to a; < 10, this problem is circumvented, 
(b) The lowest eigenstates of the stadium show some marked regularities. In our analysis, we considered 
levels only above the 50*'\ In so doing, the agreement for (p ~ Q with the GOE results (for levels up 
to 300) is very good, (c) Spectral modulations due to 'bouncing ball' orbits would make comparison 



of the model (i.e. Eq. (2.1)) with the billiard difficult. For our choice of a relatively small i/R, similar 
to the experimental situation, there are no sizeable modulations. The numerical results indicate that 
\/t — 3.87 0/00- The uncertainty in the fitting of the best Y2{x,t) is about 5% due to statistical error 
bars. Here, we determined t for a system containing approximately n = 400 electrons. The physical 
situation of interest is more likely to deal with n ~ O(IO^), and the question is whether our result will 



be changed for this region of the spectrum. 

A partial answer to this question was given by Berry and Robnik , who analysed the semiclassical 
transition parameter for a billiard threaded by an Aharonov-Bohm flux line. They found that for a fixed 
4>/(t>o, t increases as n^/^. A complete semiclassical treatment can be found in Ref. |l^, where t is expressed 
in terms of the density of states and purely classical quantities: 

Vt^ ^ . (3-5) 

where K(E) is a classical quantity which measures the average magnetic flux squared that is enclosed 
by a typical trajectory. We should like to stress that it is not the purpose of this discussion to further 
improve the semiclassical approach, but rather apply the findings of Ref. |l^ to obtain some quantitative 
results. For this purpose, we shall present the semiclassical results in a nutshell, motivate their origins, 
discuss briefly K{E) and rescale it in a convenient way. Finally, we compute K{E) for a certain geometry 
and determine the t-dependence in a way independent from the previous. 



After some manipulations on the Gutzwiller trace formula, Ref. p7| gives 1^2 (e) as 

Y,(e) = dr \r\e^^^^-^^M)/'^' (i + ^g-(AS^(.)>p... ^ J , (3.5) 

where the average (...) e should be taken over an energy interval large compared with e and the smoothing 
parameter 77 should be small compared with e. The central dynamical quantity to be computed is the 
phase space average (...)p.s. of as a function of time r, which replaces the usual average over periodic 
orbits. It is interesting to note that all the physics due to the presence of the magnetic field is contained 
in the latter quantity. Introducing the dynamical variable ip^r) = A(r) • r — eBip{T), one can write 

{AS^t))p,,, - 8(eB)2 r dt r * dt' {ip{t)ip{t + i'))p.s. . (3.7) 



/o Jo 

There is considerable work in the literature dealing with long power law tails — s- cx) of {ip{t)ip{t + i'))p s. 
and similar correlators. These works and our numerical studies assure that {ip(t)(p{t + t'))p s. decays 
sufficiently fast so that one can define 

K{E,t)= f dt'(^(0)^(i'))p.s. (3.8) 
Jo 

and (AS'^(t))p.s. w 8{eB)^K{E,T — > oo)t for sufficiently large t. The quantity K{E) has trivial E and 
A dependence, and it is convenient to single out these variables so as to obtain a correlator k defined for 
a unit area and unit velocity, which depends solely on the shape of the billiard. This can be easily done 
by a proper rescaling of the coordinates and the time; as a result, one has 



K{E,T) = \l^^n{T). (3.9) 
m 
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Inserting this into the scmiclassical expression for the transition parameter and recalling Eq. (|3 
one arrives at the relation 

^"(^) ' (3-10) 

where N^^{E) is the cumulative level density, which counts the number of single-particle states in the 
billiard up to the energy E. This result gives us the correct energy scaling of the problem, implying that 
for a fixed magnetic field the GOE GUE transition will become faster as the number of electrons n 
increases (in agreement with Ref. n9^ and our naive estimate). Moreover, it indicates that ^/i depends 



linearly on (t>/4>Q as a naive inspection of Eq. (3.1) would suggest. Eq. (3.1C) also allows one to obtain 
quantitative results by determining k for the the specific system under consideration. In Fig. 4, we present 
our result for k, obtained by a numerical average over the phase space. This was done by computing 10^ 
trajectories with randomly chosen initial points. (It is of interest to mention that such an average has 
to be performed over closed paths; otherwise the results depend on where the origin of the coordinate 
system is taken.) Our numerical result is k(t — > oo) ~ 0.06. It follows that \/t = a[N^^ (E)]^^^ (l>/4'o, with 
a ~ 1.05 in good agreement with 1.03 = 3.87/(200)^/'' obtained by the other method (where we should 
note that N'^^{E), being the number of levels, is half as small as n due to spin degeneracy). 

4 Supersymmetry Formalism 

To calculate the ensemble average g, we use Efetov's supersymmetry approach in the version described 
m Ref. 0. In that paper, the scattering problem formulated in Section 2 was worked out without the 



GOE-symmctry-breaking term in Eq. (2.1). The present calculation follows very closely the lines of 
Ref. |l|. 

The appearance of the symmetry-breaking term causes the scattering matrix not to be symmetric. 



This leads to a set of source terms not considered in Ref. 14, Moreover, the appearance of the symmetry- 
breaking term in the effective Lagrangian of the zero-dimensional non-linear sigma model derived in the 
N oo limit makes it much harder to calculate the terms of highest (eighth) order in the anticomniuting 
integration variables. These are the terms that yield the answer to our problem. (All terms of lower order 
yield vanishing contributions. This statement, well known in the absence of symmetry-breaking, applies 
also in the case considered here [^of .) There are two ways to overcome this difficulty: (i) One may use a 
parametrization of the Q-matrix developed by Altland, Efetov and lida [plf which is particularly taylored 
to time-reversal symmetry breaking, or (ii) one may use Efetov's original parametrization [p^ of the Q- 
matrix in conjunction with the method developed in Ref. ^ to simplify the form of the symmetry-breaking 
term. 
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We have tried both approaches. It turns out that method (ii) is very much simpler; The terms 
of highest order can be obtained by hand, i.e. without recourse to a computer-based algorithm. The 
reason is that in comparison with the level correlation function, the source terms needed to obtain the 
elements of the scattering matrix depend on Q in a more complicated way: The matrix Q appears in the 
denominator. In such a case, the method of Ref. |2l| becomes very unwieldy. This is why we have adopted 
method (ii). 

In what follows, we have attempted to render a reasonably complete and systematic account of the 
derivation of the explicit formal expression for g as given within the supersymmetry formalism. Though 
much of this material has appeared previously in various guises, these sources are rather fragmentary, and 
so we consider it worthwhile to present here a pedagogical treatment that is hopefully comprehensible to 
the relatively uninitiated reader who has a little background in Grassmann algebra and general random- 
matrix theory. At the same time, we aim to achieve consistency of notation and conventions. On the other 
hand, the expert reader should simply take note of our generating function as given by Eqs. (|]i|)-(|t|). 



which incorporates the new source matrix that is employed in the present calculation. This is given in 
Eq. ( [I.19| ). From here, one can proceed directly to the final integral expression over supermatrices for 
l^abp, which is presented in Eq. ( 4.42| ). 



4.1 Generating Function 

Let us consider a generating function given by 

Z(e) = T)etg-^[D + J{e)] 

= exp{-Trgln[i? + J(e)]} , (4.1) 

where the symbols 'Detg' and 'Trg' denote the graded determinant and graded trace as defined in the 
supersymmetric formalism of Ref. |l^. Here, the inverse propagator D has been extended to a AN x AN 
supermatrix, 

D = {E - H)li + ^ujL + iQL , (4.2) 

where 51 — J2c ^^"^ ^pp' ~ i~^)^^^6pp'S°'" is the diagonal supermatrix that distinguishes between ad- 
vanced and retarded parts of D. We have Dpp/ — [diag(Z), D'')]ppi where p,p' — 1, 2, with p = 1 referring 
to the retarded block and p = 2 to the advanced block. The indices a, a' = 0, 1 determine the grad- 
ing (with a ~ for the commuting (bosonic) components and a = 1 for the anticommuting (fermionic) 
components). We have also allowed for the presence of a difference uj = Ei — E2 between the energy 



arguments of the two S-matrix elements that multiply together to form g. According to Eq. (4.2), Dn 



12 



corresponds to energy Ei = E + while D22 corresponds to energy E2 = E — ^uj. We shall ultimately 
set u! 0, but keep it in the interim to facilitate comparison with other related work pl|] . 

The source supermatrix J(e) depends on a set of parameters Sm labelled by by some (multi)-index 
m, and is taken to have the general form 

J{e) = e„M„ . (4.3) 
for some set of 47V x AN supermatrices Mm(/^, i^)- It then follows that 



-Z{e) 



^ TrgM„D-i- TrgM„D-i + TYgMmD-HlnD-^ . (4.4) 

e=0 



deaden 

For the problem at hand, we shall make the choice of source matrix 



2 

= E Iimausi,W,, , (4.5) 

a,b j = l 

where 

(4-6) 

is a projector onto the p — j block with k"" ~ (— 1)"(5"" . In this case, the second term on the RHS of 



Eq. (4.4) vanishes if we differentiate with respect to £,„ — e]^^^ and £„ = e^^. Comparison with Eq. (2 
then leads to the identification 

(4.7) 



52 

\Sab\ = Q 1 Q 2 ^(^) 



e=0 



4.2 Ensemble Average 



As a Gaussian superintegral, the generating function Z{£) reads 



Z{e)^ j I?^I?^expjz^(^^,[(i? + J)^]^) 



(4.8) 



where we employ the notation 

(F(/x),G(/x))^^F(/z)G(m), (4.9) 

and V'p (/^) is a four-component supervector field. The adjoint supervector is defined by ^ = ip^ s with 

s'^P = V , = (-l)(i-")(i+P) . (4.10) 

The presence of the supermatrix s in the definition of the adjoint ensures the convergence of the final 
integral representation of the ensemble-averaged generating function as a supermatrix non-linear cr-model 
and ensures the correct combination of compact and non-compact symmetries therein p3]. 
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With the definitions 



we have the property 



where 



<y5 

SLp* 



A = 



A 

AT 



(4.11) 



(4.12) 



(4.13) 



for any supermatrix A that is diagonal in the graded indices a,a' . The eight-dimensional supervector <I> 
possesses the 'charge-conjugation' or reality property $* = (7$ which can equivalently be expressed as 
a condition of Majorana type, viz. $ — (A<l>)Tj where the matrices C, A are defined in Appendix A. In 



view of Eq. (4.12), Eq. (4.8) may be cast in a form that renders the GOE ensemble average particularly 



simple, namely 



Z{e) = / P$ expi [(D + J)*]^.) 



(4.14) 



The indices r, r' span the additional matrix structure due to the doubling of dimensionality implied by 



Eq. (4.11). We shall say that the index r labels the 'GOE blocks'. Here, we have 



D = Els, - H + \ujL + mi 



and the definitions of H and J are those implied by Eq. (4.13). In particular, 

T 



H 



H 

i/T 



(4.15) 



(4.16) 



where the matrix is defined in Appendix A. We should note the following symmetry properties: (i) 
Transformations T : $ i-^ r$ that preserve the relation $* — C$ obey 



T* = CTC- 



and (ii) transformations T : ^ ^ T$ that leave the bilinear invariant are those which satisfy 



(4.17) 



(4.18) 



The set of 8x8 supcrmatrices respecting conditions (4.17) and (4.18) forms a supergroup with compact and 
non-compact bosonic subgroups. In the absence of symmetry-breaking and source terms, the generating 



function in Eq. (4.14) is invariant under this supergroup. 



Special attention should be paid to the induced r-space structure of the source supermatrix J: 

2 2 



a.b j=l a,b j=l 



(4.19) 
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where e^^''"' and £)j^'''' are the parts of e^j^ symmetric and antisymmetric in the indices a, b respectively, and 
we have introduced = <E) I2, l2(i) = ® t^- Consequently, the symmetrized form of \Sab\'^ 
can be written as the sum of two terms: 



Q2 



Q2 



e=0 



4 r,{a)\r,{a)2 



Z(0, 



(4.20) 



£=0 



^'^ab ^^ba 

for a ^ b. Cast in this form, the contribution to |S'ahP from the first term on the RHS of the equation 
above is generated by the essentially same source matrix as used in Ref. |lj for the pure GOE problem; 
the second term is generated by a new source matrix. This latter term must vanish for i = and accounts 
for the asymmetry of the S-matrix. 



The factor in the integral of Eq. (4.14) which undergoes ensemble averaging can be decomposed as 

■■ H « A 

(4.21) 



exp{-^ ^($p„(H$)^,,)} ^exp{-zTri7(^)5(o)} • exp{ v/tT^Tr (^'5(3) } 

p,a,r ^ ^ 

where we have introduced the ordinary N x N matrices in the level indices /x, i^, 



and we take = 1. Gaussian ensemble averaging leads to the relations 

-s r 



exp|-iTri/(^)5(o)} expj-— trg5' 

— f — ; — { f ^ 

exp|yt7^Tri7(^)5(3)| = exp <^ - — — trgr^S'r^S' ^ 



where 



is a supermatrix. These results are easily derived with the aid of the relations 

^5(«)(.,a.)5(°)(m,^) = itrg5^ 
^5(«)(z.,m)5(")(^,A.) = itrg5^ 

fi.y 

and 

^5(3)(^,^)5(3)(^,^) = +ltrgr'Sr'S, 
^ 5(3) (z.,^)5(3) (,.,;,) = -itrgr^^r^^. 
The ensemble-averaged generating function then becomes 

Z(i) = j I?$e'^*(^) expj^ ^($p,,[(^l + ic^i + ir!L + J)$];^, 



p^a,r 



where 



iC*{S) = - 



trg + —tvgT^ST^S 
N 



(4.22) 



(4.23) 



(4.24) 



(4.25) 



(4.26) 



(4.27) 



(4.28) 
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4.3 Hubbard- Stratonovich Transformation 



Following Ref. ^ let us introduce unity under the ^-integral in the form 

1 = J I?(Te*^("'^) , U{S) = J Va^^^"-^^ , (4.29) 

followed with an interchange of the cr and $ integrations, where a spans some appropriate space of 
supermatrices, not yet specified. The quantity W{a, S) can be any function provided the integral in 



Eq. (4.29) is convergent. The choice 



W{a,S) = -C* [S-'-^a 



(4.30) 



leads to 



i£*{S) + iW{a,S) 



N 



trgcr 



N 



trg T^CTT^a 



y ^($;,[(S + |^r3Er3)$]P, (4.31) 



where S^j^ = aS^^. The quartic dependence on <i> in the exponent of the generating function implied in 
£*(5) is thus eliminated, and if the space of cr's is chosen such that M{S) = JV{0) (i.e. shifts are allowed), 
then the remaining quadratic dependence on <i> is amenable to exact integration, which yields 

t 



Z{e) 



Va exp<; - — 



trgcr^ + — trgr'^crr^CT 



t 



exp I -\ Trg In E\ + \uL + i9.L + J - A(S + — r^Er^^ 



The appropriate space of matrices a is given by those which assume the form 



a = T-^{P -irjL)T 



(4.32) 



(4.33) 



where the P, which represent massive modes, are essentially Hermitian and block-diagonal in the p- 



index, and the T span the set of matrices defined by Eqs. ( 4.17 ) and Eq. ( 4.18 ) modulo the subgroup 
that commutes with L. The positive constant 77 is required for convergence, and its precise value will be 
determined from the saddle-point equation. Next, we write P — irjL = (Tq + 5P , where <jq is the unique 
diagonal saddle-point that lies within the integration manifold, and 5P are the fluctuations around it. 
Thus we have a = <tg + T^^SPT with ctg T^^a^T . The massive fluctuations can be integrated out 



in the limit N 00, which leads to an expression for Z{e) identical with Eq. (4.32) except for the 
replacement cr ^ ctq everywhere. 



With E ~{) and neglecting u) (which is assumed to be of order 0(A^^^)), the saddle-point equation 
reads ctq + Cq"^ = 0. Setting = x\ + yL leads to a; = 0, y = ±i. According to Eq. (4.33), the correct 
choice is ctq = —iL. So we write ctg = —iQ where Q = T^^LT. Having set E = 0, we can write 



Z{e) = JvQe 



i£otf(Q)gj£=ouroc(Q;J) 



(4.34) 



16 



where 



and 



with 



«Aourco((3; J) = Trgln 



itji + iVlL + i AQ 



(4.35) 



(4.36) 



(4.37) 



The summation over the level index /i (implied in the trace 'Trg') can be performed by virtue of the 
orthogonality relation of Eqs. ( 2.10| , 2.1l|) , from which one can derive the result 



Tr(w + = {N- M)oj'' + ^{Ld + X^Y , 

II. ' ^ 



(4.38) 



for n = 1, 2, . . . and any constant lo. (In our case, w = ~iuj /2.) Consequently, Ccs{Q) can be decomposed 
according to £off(Q) = ^hcc{Q) + 'Cch(Q) with 



i£frcc(Q) 



- trg r^Qr^g- y trgln(icjL + iAQ) , 



1 



1 



(4.39) 



icjL + iXQ 

Next, we expand C^cdQ) and £ch(Q) in the small quantities uj and i/A^ to leading order. Using the fact 
that — 1, this yields 

iC,^{Q) = -i^trgln(^l + ^LQ^ . (4.40) 

We note that iC^cdQ) coincides with the analogous quantity derived in Ref. ^ At this stage, we 
implement the assumption of equivalent channels with ideal couplings, as already discussed in Section 2, 
by setting Xc = A for all c. Thus we obtain 



M 

iC,i,{Q)^^—iYg\n{l + LQ) , 
noting that 1 = Af by definition. The source term can be treated similarly. 



(4.41) 



Letting w — > 0, we finally arrive at the expression 



\Sab? 



-y trgln(l + QL)~^- trg gr^Qr^ 



(4.42) 



VQ exp 

2 

tro- , ^ , 

QL '1 + QL 

The space of the supermatrices Q can be parametrized by 16 real variables. These can be chosen to 
comprise 2 'azimuthaP angles, 3 'Cayley-Klein' parameters, 8 anticomniuting variables and 3 eigenvalues 



.^trg 



1 1 (1) 1 1 (2) 
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(1 compact and 2 non-compact). All but the three eigenvalue parameters can be integrated out analyt- 



ically. The final threefold integral giving the exact result for |S'afcp in its maximally reduced analytical 
form must be further evaluated numerically. Our result in Eq. (4.42) differs from that of Ref. |l^ by the 
appearance of the new source term containing l2(p) and the symmetry-breaking term involving in the 
exponent. 



5 Exact Results 



Using the method developed in Ref. ||, we can perform in Eq. ( 4.42 ) the integration over all the anti- 
commuting as well as over several of the commuting variables used in the parametrization of Q. We are 
then left with a threefold integration over real variables (the 'eigenvalues' in Efetov's parametrization). 
A summary of the technical details involved in the derivation is presented in Appendix B. 



For ease of subsequent discussion, it is convenient to separate the exact result for jS'abp into two 
contributions: 

=\S^ +\S^ , (5.1) 



where |S'ip represents the part arising from the first source term (j = 1) of Eq. (4.42) and |S'2p that 
from the second source term (j = 2). Then we have, for a 7^ 6, 



where 



And, 



d\ / dXi 



1 + A 
Ai + A2 



M 



1-A \l + Xl-2 



1 + A (Ai + A2)2 



1 

At 

n 



\l - A2)(1 + e-2*(i-^')) - (A? - Xl){l - e-2*(i-^')) 
(1 - A^)e-2*(i-^') + e-2*(i-^'))" 

n = x^ + + 2AA1A2 - 1 



1 



1 + A 



-fl nOQ nOQ 

2./_, "V, "■^'VA. + A, 



M 



-2t(\l-l)_ 



n 



1 „ g-2i(l-A") 



1 - A 



i(l-A2) y (l + A)(Ai+A2) 
(Ai - A2) [3 + |(1 + A)(AiA2 - A)] + 4t [(1 + A)(A2 - AAi) - Xl{Xi + A2)] 



+ 



(l_e-2t(i-A^) 
At 



l + -(l-A)(AiA2-A) 



Ai + A2 



[(l-A)(A2-AAi) + A2(Ai-A2)] 



(5.2) 



(5.3) 



(5.4) 



18 



We remark here that (i) for t = 0, |<S'2p = 0, so that only jiSiP contributes to the GOE Umit, while (ii) 
for t — > oo, IS'ip = |'S'2p, i-e. the two terms gives equal contributions to the GUE limit. 



The S-matrix result above can be converted into a conductance by virtue of Eq. (1.2). We define the 
quantity Sg to represent the positive deviation of the mean conductance g from the GUE limit t oo, 



6g{t, M) = g{t ^ M) - g{t, M) 



(5.5) 



5.1 Analytical Results for Special Cases 



For general values of t, one must resort to a numerical evaluation of Eqs. (5.1)-(5.4). However, various 
special cases can be treated analytically. Here, we shall quote the results, relegating details of the 
derivations to Appendices C and D. (a) In the GUE limit {t oo), one integration drops out to yield 



1 



(Ai-A)= 



1 - A Ai - 1 



1 + A Ai + 1 



The corresponding value of the mean conductance in this limit is then 

9gve[M) = Y , 



(5.6) 



(5.7) 



in agreement with Eq. (U8). The integral in Eq. ( p.6[ ) does indeed correctly represent the GUE limit 



of I .Sab p. This was worked out separately along the lines of Ref. 14, replacing the GOE Hamiltonian 
ensemble used there by a GUE Hamiltonian ensemble, (b) In the GOE limit {t — Q), one has, for a ^ b, 



cl\ 



M 



1- A2 



1 + A 

Ai + A2 ; (A2 + + 2AAi A2 - If 



1-A \l + \l-2 
1 + A ^~ 
1 



(Al+A2)2 



71/ + 1 

This leads to a mean conductance for zero field that is given by 

M2 



5goe(M) 



(5.8) 



(5.9) 



2(M+ 1) ■ 

We have also checked that, after a suitable change of integration variables, the integral formula above 
agrees with Eq. (8.10) of Ref. m upon setting Tc — 1 for all c and considering a ^ b. Details are provided 



in Appendix C. (c) One can also perform a large-i expansion of Eqs. (5.1)-(5.4). For values M > 2, one 
obtains 

= i - T7?-^7^ + O f . (5.10) 
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The corresponding result for the large-i tail of the weak-localization term is 



5g{t,M) 



M3 1 



O 



M2 - 1 16< 

The derivation of this result, which is not entirely straightforward, is outlined in Appendix D. 



(5.11) 



5.2 General Magnetic Field 



For general values < t < oo, the eigenvalue integrations in Eqs. ( |S.2| , 5.4) can no longer be carried out 
analytically. However, progress can made examining the implications of assuming that (5.g(i, M) is indeed 
Lorentzian in the magnetic field: B ^ 4'/4'o ^ Vt- Accordingly, we write 

1 _ 1/a 



6g{t,M) = 



(5.12) 



a + bt 1 + {b/a)t ' 

The unknown parameters a and b can be determined completely from the three special cases discussed 



above. To find a, we set t = in Eq. (5.12). This gives 

- = <5,g(0,M) 



g{t oo, M) - g{0, M) 
M 



2{M + 1) 



(5.13) 



with the help of Eqs. (p]q) and (5.6). To find b, we expand Eq. (5.12) in powers of 1/t, which yields 



5g{t,M) = 



bt 



1 



l + a/{bt) 



1 

bi 



(5.14) 



Thus 1/6 is the coefficient of 1/t in 5g{t, M) as given in Eq. ( 5.11 ). Consequently, we see that 

b 8(M-1) 



a M2 

Substitution of Eqs. ( |5.13| ) and ( |5.15| ) into Eq. (|1|) gives us the result for the expected Lorentzian 

M 1 



(5.15) 



6g{t,M)^ 



2(M+1)^ ^ 8(M-1) ^ 



(5.16) 



M2 



6 Numerical Analysis 



In order to verify whether Eq. ( |3.16 ) is indeed a good representation of Eqs. (5.1)-(5.4), we should 



compute jiSabp numerically as a function of t and M. This data is most easily digested if we plot the 
curves of reciprocal weak-localiztion versus t, i.e. 



vs t , 



(6.1) 
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for various values of M. According to the Lorcntzian hypothesis, we expect to find 

5g-^{t, M) = a{M) + b{M)t , (6.2) 
i.e. a straight hne in t for each value of M . 



Unfortunately, the exact integral expression for \Sab\'^ as given by Eqs. ( |).lD -( [5^ ) is not directly 
amenable to numerical evaluation. It must first be subjected to a number of changes of integration 
variable in order to improve the behaviour of the integrands. This procedure is outlined in Appendix E. 



The numerical results for the curves specified by Eq. (|6T|) are presented in Figures 5 and 6. One can 
make the following observations: (i) For M = 2, 3, 4, 5, . . . and t ^ 100, the numerical plots agree with 
the 'theoretical' straight lines to better than 1 part in 10^, and the degree of agreement improves as M 
is increased, (ii) Odd values of M are unphysical, though the threefold integral remains well-defined. 
This holds true even for M — 1. We observe, however, that the M = \ curve deviates significantly from a 



straight line. Moreover, Eq. (5.16) does not even make a sensible prediction for what straight line could 
have been expected in this case. To illustrate the extent of the deviation, we plot in Figure 7 the weak- 
localization term 5g as a function of \/i for i\/ = 1, and compare it with the corresponding Lorentzians 
obtained from fitting to the large--\/t tail of the M — 1 curve, and also to its FWHM. The inference 
that one draws from all this is that, viewing Sg{t, M) as a continuous function of the real variable M, a 
transition to almost Lorentzian behaviour happens somewhere in the region 1 < AI < 2. 

One technical point, especially relevant to the case M = 1, that is worth mentioning is the following: 



For given M, the quantity \Si\^ varies monotonically from the GOE value of 1/(A/ + 1) at t = to half 
the GUE value, viz. 1/(2A/), in the limit t oo. Thus for M — 1, it should remain constant with a value 
of i in spite of its complicated structure in the variable t. This is precisely what is observed numerically. 

7 Discussion 

To determine the magnetic-field scale on which the GOE GUE transition occurs, let us look at the 
full width at half maximum (FWHM), which we shall denote by Fi/2{M), of the 6g versus \/t curve: 



6g — l/[a + 6(Vt)^]. From Eq. ( 5.16 ), it is easy to see that 

M 
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Thus we conclude that the GOE GUE transition is delayed as the number of open channels increases. 



We should also observe from Eq. (5.16) that the limits t — > oo and Af — > oo do not commute: 



\im Sg{t,M) = for all Af , 



lim 5g{t, M) = - for all t . (7.2) 

Af^oo 2 

In the former case, one trivially remains in the GUE limit whatever the value of M ^ while in the latter, 
no crossover is observed, with the GOE weak-localization term holding for all values of t. 

This phenomenon can be related to two extreme cases in problems involving quasi-one-dimensional 
rings. When one considers an isolated quasi-one-dimensional ring (hence not coupled to external channels) 
such as in the calculation of persistent currents , the transition from GOE to GUE sets in very rapidly 
with increasing flux angle — at a value around 4>/(j)Q ^ [d/ {ttEc)]^/'^ , where Ec denotes the Thouless 
energy. On the other hand, in the problem of AAS oscillations in the same quasi-onc-dimensional 
ring now coupled to two leads carrying a large number of channels M , one computes the weak-localization 
term as an expansion in terms of diffusons and Cooperons Here, one does not observe a transition 
to GUE behaviour for any value of the flux angle. The oscillations of the magneto-conductance persist 
for all values of (j)/ (po, as the Cooperon (a periodic function in 0/0o) is never suppressed. The hypothesis 
that the GOE GUE transtion is delayed with increasing channel number clearly renders these two 
distinct observations consistent. 



The broadening of the Lorcntzian with increasing channel number implied by Eq. ( ^.16 ) can be 



simply understood in terms of the two relevant time scales that enter the problem, (i) We have r^ix — 
the average time taken by the time-reversal violating perturbation (due to B) to mix the GOE levels. 
The width associated with this time scale is 

" t 



2tt 

~d [N 

2'Kt A2 



rr(A)|2 



in which case 



(7.4) 



r ■ 2fd 

^ mix ^lyUj 

since d — ttX/N. Thus, a stronger magnetic field requires a shorter time to mix the levels, (ii) We also 
have T(jcc — the decay time of states in the internal region due to the coupling with the leads. The 
associated width is given by 

dec = ^ , (7-5) 



and this leads to 



Tdcc = f; — = -rrr > (7-6) 
Idee Aid 
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from which we see that for more open channels, the electrons spend less time within the stadimii. There- 
fore, an increase in the number of open channels implies a decrease in Tdoc, which means that the same 
degree of time- reversal violation will require a smaller r^ix', and this is achieved by a stronger magnetic 
field. The ratio rmix/rdcc ^ t/M provides a measure of the time-reversal violation, and we note that for 
M 3> 1, 5g depends on t only via the combination t/M. This feature is not explicit in the semiclassical 
result of Ref. |l^ but can be inferred indirectly ^ . 

If one considers the problem of calculating 5g{t^ M) in the framework of the asymptotic expansion 
of Rcfs. 1^ and|^, which essentially proceeds in inverse powers of M, then one is able to identify analogues 



of the diffuson and Cooperon as Hp = M and 11^ — M + 8t, respectively. Eq. (5.16) has a very simple 
interpretation in terms of these quantities: If one considers M to be large (compared with unity) but 
t/M not necessarily small, then one can write 

6g{t, M) 2 (^) (M.) none , (7.7) 



M»i \ 2 y V 2 , 

where the first three factors represent the spin degeneracy and the summations over the incoming and 
outgoing channels. In this approach, a diffuson factor is always implied by the structure of the lowest- 
order contribution to the source term, and reflects the coupling with the leads. With increasing t, the 
Cooperon contribution is damped out, and with it Sg. 

We have treated the stadium billiard strictly as a one-particle system in which classical chaos and 
magnetic flux combine to cause the GOE GUE crossover transition. In actual fact, we should have 
treated the stadium as a problem involving a large number (say A) of non- interacting electrons. In such 
a system, the eigenvalues of the many-body Hamiltonian H — J2f=i ^U) '^ill '^'^^ obey GOE (or GUE) 
statistics, even if the eigenvalues of each of the single-particle Hamiltonians H{j), j = 1, . . . , A do. Does 
this fact necessitate a rethinking of our method and of our results? The answer is no. Indeed, using the 
argument of Ref. one can easily show that for non-interacting fermions, it is legitimate to calculate 
the single-particle scattering matrix without paying attention to the presence of the other fermions. 

Finally, let us attempt a comparison of our theoretical result ( ^.16| ) with experimental data. There 



are several reasons why a direct comparison of our Eq. (5.16) with the results of Ref. ^ for g is not 
possible. First, the stadium used in the experiment possesses a twofold reflection symmetry. Hence, in 
the closed stadium there exist classes of states pertaining to different quantum numbers with respect to 
this symmetry. This fact is not taken into account in our modelling of the Hamiltonian. Second, the 
stadium is known to have a family of 'bouncing ball' orbits which might obscure the connection with 
our random-matrix model. In view of the way the leads are attached to the stadium, it is likely that 
these orbits affect the experiment only slightly. Thirdly, the data of Ref. were obtained by averaging 
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over the gate voltage or, equivalently, over the number of channels and the area of the stadium. Finally, 
the data show a marked dependence on temperature which we do not take into account. For purposes 
of orientation, we nonetheless take the data of Ref. ^ at face value. For large B, the value of g is close 



to 3. This suggests M « 6, recalling that 5gue(-^^) ~ M/2. Insertion of this value of M into Eq. (5.9) 
for the zero-field value of 5g yields a theoretical prediction of 0.43, to be compared with a value of 1/3 
experimentally. Also, one can write 

Vt = k(j)/(l}o = {0.2A)kab , (7.8) 
where a is the area of the stadium measured in (/im)^ and b is the magnetic field in units of mT. With 



the aid of Eq. ( |7.lD , this gives the theoretical FWHM in 6 to be 



1 / M2 

for M = 6, a ~ 0.5 (/im)^, and taking k w 5.8. This value of k results from estimating the number n of 
electrons in the stadium to be about two thousand and using the scaling relation 

fc4oo = 3.87(-) , (7.10) 

as discussed in Section 3. Both the values for the height of the Lorentzian and the FWHM differ somewhat 
from what we read off the experimentally fitted Lorentzian shown as a dotted line in Figure 3(b) of Ref. 0, 
namely 1/3 and 5 mT, respectively. The discrepancy is consistent with the presence in the experimental 
data of effects, such as finite temperature and inelastic scattering, that will tend to wash out coherent 
phenomena, and manifest themselves as a fiattcning of the peak and broadening of the width. 



8 Conclusions 



We have derived a universal expression for the flux and channel-number dependence of weak localiza- 
tion, Sg, in classically chaotic ballistic microstructurcs. For AI > 2, 6g{t, M) is almost Lorentzian in 
\/t ^ B. The FWHM of 5g{t, M) versus \/t grows with channel number essentially as ■y/Af/2, so that the 
GOE GUE crossover occurs at t ^ A//8. This fact is arguably the most interesting result of this 
paper: An increase in the number of channels causes an increasing delay in the GOE — > GUE crossover 
transition. 

It is remarkable that for values of M as small as two, 5g{t, M) is almost Lorentzian in iJ. This 

is reminiscent of results for the S- matrix auto-correlation function versus energy in the GOE limit. In 
Ref. it was found that this function is close to Lorentzian even for small M provided all channels 
have Tc — 1. For 'inequivalent' channels, i.e. channels with different Tc's, marked deviations from the 
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Lorentzian form occur for M as large as ten ||lO[|. We speculate that, similarly, gating one of the leads 
connected to the stadium may significantly affect the functional dependence of 5g{t^ M) on t. 
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Appendix A 

In this appendix, we collect the various constant matrices that appear throughout the paper. Matri- 
ces arising from the graded structure of the theory and the symmetry-breaking due to the presence of 
advanced and retarded components are given by 

fc"" = ( -1 ) ' ^pp' = ( -1 ) ' " ( -fc ) • 

V / aa' V / pp' \ / pp' 

Matrices associated with the doubling of dimensionality in the GOE problem relative to the GUE are 
given by 

A,.' =(l \ \ , C^s(g>A. (A.2) 

\ / rr' 

Matrices arising from the presence of magnetic field and the source terms are given by 

Trr' = ( -1 ) ' ^PP'W ^ ( ) ' ^PP'C^) ^ ( -fc ) ■ i'^-'^) 

\ / rr' V / pp' \ / pp' 

In all cases, the explicit indices indicate the space to which the displayed block structure pertains. 

Appendix B 

In Efetov's parametrization, one writes Q = UDU, where 

/ cos ^ i sin ^ \ n f (^B \ f-r, ,^ 

^^[~^sin9 ~cos0 j ' 0^ ) ' ^^'^^ 

with 
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and we have the relations u = ,v — kv^k. The rows and columns of matrices of dimensions 8, 4 and 2 
are labelled by the indices (p, a, r), (a, r), and r, respectively; the indices follow in lexicographical order. 
(With a = 0, 1, this order differs from that of Efetov The 4x4 matrix u is itself expressed as a 

product of two matrices u ~ uiU2 (whence u — U2U1) where 

/ 1-27777 + 6(7777)2 -2(1-27777)77 \ „ - f Fi \ , . 

"i^i 277(1-27777) 1-27777 + 6(7777)2 ] ' "2-1 F2 J ■ ^"■■^) 



Here, 



and 



"^l "^i ] , V = V^= "^l ^2 ) , (B.4) 



Fi = exp(i(/.T3) , F2^i^_^* J j , (B.5) 

where jwp + |z|2 = 1. We note that F2 is an SU(2) matrix, which can be expressed in terms of uncon- 
strained parameters according to 

- = 'f"^: „. (B.6) 

1 + 77!,^ + |77T.l 1^ 1 + 777^ + |771l|^ 

Similarly, we write v — V1V2 (whence v ~ V2V1) where 

_/ l + 2pp + 6(pp)2 -27(l + 2p/7)p \ „ \ , . 

27p(l + 2pp) l + 2pp + 6(pp)2 ) ' "2-1 $2)- ^^■^) 



Here 



and (f>i = exp(7X''" )) ^2 = 1- The normalization of it's and t;'s is given by uu = vv = 1, UjUj = VjVj = 1, 
for j = 1,2. 

In the parametrization given above, Oi, 62, 0, m, (j), x are real commuting variables, ttt-i is a complex 
commuting variable, while 771, 77J, 772, 772, pi, Px, p2, P2 are all anticommuting variables. If we also introduce 
Ai — cosh 01, A2 = cosh 02 and A — cos0, then in terms of these final variables, the integration measure 
on the space of Q-matrices is given by 



2(1 — A2) 

— /r.. dXidX2dXVF2 d(l)dxdr]idrildr]2dr]2dpidpldp2dp2 , (B.9) 
(647r/c)^ 



where TZ = Xf + X2 + X — 2AA1A2 — 1 and 'DF2 denotes the SU(2) integration measure 

4 

'DF2 = 5—^ -r-^d771 d(Re 7711 ) d(Im 7771 ) . (B.IO) 

7r^(l + 777^ + iTnip)-^ 



One should note that Eq. (B.£) assumes Efetov's convention ||T^ for integrals over the anticommuting 
variables, / dxx — ^1 rather than the convention of Ref. |lj, namely J dxx = (27r)^^/2 xhe integration 
region corresponds to 1 < Ai, A2 < 00, — 1 < A < 1, —00 < 777,Re77ii,Im77ii < 00, < (p^x < 
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To perform the Grassmann integration, we recall that only that part of the integrand which is of 
the highest (eighth) order in the anticomniuting variables contributes . The Af-dependent term of the 
integrand depends solely on the generalized eigenvalues 9i,92,9. By diagonalizing the matrix 9, we find 
that 

exp {-iAf trgln(l + QL)}^ Detg-*^/2(l + QL) 
After making the substitution 



1 + cos 6* 
cosh 9i + cosh 6*2 



M 



(B.ll) 



(l + QL)- 



1 zsin6'(l + cos6')"i 

isin^(l + cos^)"^ 1 



U 



(B.12) 



we find, for the symmetry-breaking and source terms, the explicit forms 
trg(Qr^QT^) = trg[(cos ^ U2-wit^wi-«2)^ 

+ 2 sin ^ U2-UlT^Ui-U2 Sm9v2-VlT^Vi-V2 + {cos 9v2-ViT^ Vi-V2)'^] , 

trg[(l + QL)-i/i(l)(l + QL)-'h{2)] = 

— J trg[(l + cos^)^"'^ sin^U2-uiA;ui-U2(l + cos^)""'^ sm9v2-vikvi-V2], 
trg[(l + QL)-i/2(l)(l + QL)-i/2(2)] = 

— J trg[(l + cos^)^"^ sin^u2-uiA:T^ui-U2(l + cos^)^"^ sin ^ W2-wiA;T3t;i-W2] 



(B.13) 



(B.14) 



(B.15) 



The anticommuting variables enter the traces in Eqs. ( B.13D -(B.15) only via the matrices ui^ui and 



vifivi, with /i = T^,k,kT^. By employing Eqs. (B.3)-(B.8), we find that these matrices have the form 

uifiui = UffifUj , vifivi = Ug^igUl , (B.16) 
where Uf, Ug are the unitary matrices {UfUj = UgU^ = 1) defined by 

Uf=\n X 1 , [/, 

with 



1 

exp(7/) 



1 

exp(7g) 



(B.17) 



7/ 



-riiT]2 



V1V2 







7s 



P1P2 

-pIp*2 



(B.I8) 



We note that (i) the matrix Ug commutes with 9, V2 and V2, (ii) the matrix Ty depends solely on 1^2, r]2, 



and (iii) the matrix depends solely on P2,P2- Making use of Eq. ( B.16| ) and (i), we rewrite the traces 



standing on the RHS of Eqs. (B.1S)-(B.15) as follows: 



tYgiQT^Qr^) = tYg[{cos 9 Ulu2UfTfu}u2Ugf 



(B.19) 



+ 2sin^C/tu2C//T)C/]u2C/gSin^iJ2Tg3u2 + {cos 9 V2T^ V2)'^] , 
trg[(l + QL)-i/i(l)(l + QL)-'h{2)] = 

-jtTg[{l + cos9y'^ sin9Ulu2UfkfUju2Ug{l + cos 9y^ sin 9 V2kgV2] , 
trg[(l + QL)-i/2(l)(l + QL)-i/2(2)] = 

-i trg[(l + cos9)-^ sm9Ulu2UfikT3)fUju2Ug{l + cos9)-^ sm9v2ikT3)gV2] ■ (B.21) 



(B.20) 
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The calculation can be greatly simplified by using a trick introduced in Section 4.1 of Ref. |5|, which 
involves passing from the SU(2) commuting variables specifying U2 to new (primed) SU(2) commuting 
variables specifying u'2 defined by 

u'2 = U]u2Ug , u'2 = Ulu2Uf , (B.22) 

with u'2 being the same function of the variables (f),m',m'i as U2 is of the variables TOjTOi. The new 
commuting variables m', m'^ are functions of the old commuting variables m, mi and of the products ?7i?72, 
r]l'i]2, P1P2 and plP2- The Berezinian of this transformation of variables is equal to unity. For simplicity, 
the primes will be suppressed in the sequel. We note also that the traces depend on the anticommuting 
variables only via the matrices fi / and fig. Furthermore, the exponential symmetry breaking term depends 
only on the anticommuting variables ?72 7 ?72i P2, P2- 

trg((5r^Qr^) = trg[(cos ^It2T)u2)^ + 2 sin^U2r)u2 sin0 W2r^t;2 + (cos ^TJ2T^W2)^] 

= 8(sinh2 02 + \z\^ sin^ 9) + 64:{rj2V2 - P2P2)(cosh2 62 - cosh^i cosh 6*2 cos 6* - |zp sin^ Q) 
— 64[772P2 cosh 6*1 sinh(?2 sin6'z* exp(— i(</) + x)) + c.c] 
+ 64[772/02 sinh6'i cosh02 sin0 w* exp(— i((/i — x)) + c.c] 

+ 256?72?72P2P2[cosh2 ^\ - cosh^ 6*2 + (|wp - \z\^) sin^ d\ . (B.23) 

Thus, 

exp{-inrg(QT3QT3)} = exp{~2t{smh^ 62 + Izpsin^g)} 

•{1 — 16t{r]2r]2 — /92/52)(cosh^ O2 — cosh 01 cosh 6*2 cos0 — |zp sin'^ 0) 
+ lGt[ri2P2 cosh 01 sinh02 sin 9 z* exp(— i(0 + x)) + c.c] 

— l&t[rj2P2 sinh^i cosh 02 sa\9w* exp(— i(0 — x)) + cc] 

- 64:tT]2V2P2P*2 [cosh^ 01 - cosh^ 02 + (|wp - |z|2) siu^ 
+ 4i(cosh^ 02 - cosh 01 cosh 02 cos0 - |zp sin^ 0)^ 

+ it sin^ 0(|u;|2 sinh^ 0i cosh^ 02 + \z\^ cosh^ 0i sinh^ 02)]} . (B.24) 



For working out the integral, only those parts of the source terms which, on multiplying Eq. ( B.24 ), 
yield terms of the eighth order in the anticommuting variables are of interest. The explicit forms of these 
parts are as follows (cp = contributing part): 

trg[(l + OL)-i/i(l)(l + QL)-i/i(2)] 

= — i trg[(l + cos0)~^ sin0?l2fc/M2(l + cos0)~^ sin0TJ2/cgW2] 
^ -327yi?7*pipJ(l - 8??2r/2)(l + 8P2P2) 
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•[(cosh^ 6*1 + cosh^ 62 - 2)(cosh6li + cosh 6(2)"^ + sin^ 6 [l + cose)-^] , (B.25) 
trg[(l + QL)-^h{l){l + Qiy^hm 

= -itrg[(l + cos^)"^sin^?l2(fcr^)/W2(l + cos^)"^sin0U2(fcr^)gU2] 
-32?7ir;5'pip^(cosh6ii + cosh6l2)"^(l + cose)"^ 

•{(1 - \&mV*2){^ + 16p2P2)(cosh6li + cosh6'2)(l - cos6l)(|u'p - Izp) 

+ (cosh 01 — cosh02)(l + cos^^) + [I6772P2 sinh02 sin^^ z* exp(— 1(0 + x)) + c.c] 

— [I6772/O2 sinh6'i sinf? w* exp(— i(0 — x)) + C-C.]} . (B.26) 



Using Eqs. (B.24)-(B.26), the desired leading part of the integrand (i.e. eighth order in the anticomniuting 
variables) can be found. This part is independent of <j) and x, and its dependence on the SU(2) variables 
is carried by the factors x"exp(aa:), where a — tsir? 6 and x — \w\'^ — (The index n ranges over 

n = 0, 1, 2 for the source term of Eq. ( |B.25 ), and over n = 0, 1, 2, 3 for the source term of Eq. ( B.26| ).) 
The integrals over the anticommuting variables and over the angles (j) and x can now be easily done. In 
particular, the integration over the SU(2) parameters can be carried out with the help of the formula 
(Section 4.2 of Ref. |) 

r (7" / sinln n \ 

(B.27) 



da" \ a 



for n = 0, 1, 2, 3. Having performed these integrations, we finally obtain \Sab\^ hi the form of Eqs. (5.1)- 

(ph. 



Appendix C 
GUE Limit 



We begin our discussion of the GUE limit (t — > 00) by outlining the derivation of the integral appearing 



in Eq. (5.(:). For l^ip, the contribution that survives the large-t limit in Eq. (5.2) comprises the term 
whose coefficient is 4i, but neglecting the exponentials e~^*(^~'^ \ i.e. 



1^1 1 



t—^OO 2 



-1 /'OO 

dX 



1 + A 
Ai + X-2 



M 



-2i(A^-l) 



1 - A Xl + Xl 



1 + A (Ai + A2)2 



4a| 
n 



For \, this integral is dominated by the region A2 ~ 1. Thus we write A2 = 1 + z, so that 



(C.l) 



(C.2) 



The rest of the integrand should also be expanded in powers of a;, retaining only the lowest order contri- 
butions. This is equivalent to the replacement A2 — + 1 everywhere except in the exponential on the RHS 
of Eq. (|a]). We note that 

TZ - (Ai - Xf . (C.3) 
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The exponential integral over z can now be trivially carried out, and we obtain exactly half the integral 



expression of Eq. (5.6). The contribution to |S'2| in Eq. (5.4) that survives the large-i limit consists of 
the terms with coefficient 4t, upon neglecting the exponentials e^^*^^^'^^). This can be simplified to 



+1 



d\ \ dXi I d\2 I ^ I e 



Ai+Aay 7^ (l + A)(Ai+A2) ■ ^ ' 



To extract the large-t limit completely, the procedure used above for can be employed; namely, 

we write A2 = 1 + z and expand in z except for e~'^*^. This also yields half the integral expression of 



Eq. (5.6), and hence the result is established. 



To now evaluate the GUE integral of Eq. (^^) analytically, we first cast it in the equivalent form 



I^gueP - 2 / d\ d\i 



-1 



l + Aiy Ai - A (1 + A)(l + Ai) 



Then we make the change of variables 



A 

Ai 



1 — w 

1 + w ' 

1 — w 2x 



1+w 1+w 



The corresponding Jacobian is given by 



and we also note that 



9(A,Ai) ^ 4 
d{x, w) (1 + w)-^ ' 



1 + w 



Ai - A 



1 



2x ' 

_ {1+wf 

(1 + A)(l + Ai) " 4(1 + 2;) ' 
1 + A 1 



1 + Ai 



1 + a; 



Putting all this together, we obtain 



5GUE| 



dx / dw — 

Jo x{l + xY'+^ 



1 

M 



(C.5) 



(C.6) 



(C.7) 



(C.8) 



(C.9) 



GOE Limit 



To bring the integral in Eq. (p.8[) into the form that appears in Ref. h4l we apply the transformation 



l^2^l, 



Ai 



(1 + ^l)(l + ^2) + + 2^^1(1 + ^1)^*2(1 + M2) 



(1 + ^i)(l + /i2) + MiAi2 - 2a/mi(1 + Mi)a^2(1 + M2) 



1/2 



1/2 



(C.IO) 
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after re-expressing the non-compact integration in Eq. ( JS.q ) as one over the triangular region A2 < Ai < 00, 
1 < A2 < 00, so that Ai — A2 > everywhere. The associated Jacobian is given by 

i9(A, Ai,A2) ^J.l-^J.2 



and we note that 



This allows us to write 



i9(m,Mi,M2) + Mi)M2(1 + Ai2) ' 

7^ = 4(/x + /il)(/l + /i2) . 



(C.ll) 



(C.12) 



M 



- f dn f dui f duo — rnrr 

8 Jo Jo Jo v/MTT7IIWTT^[(i + Mi)(i + ^2)] ^ / 



(M + Mi)^(M + Ai2)2 



1 



1 



1 - /.i 1 + /ii 1 + ^J.2 



(C.13) 



which agrees with Eq. (8.10) of Ref. |Tj for a 7^ b, having set Tc — 1 for all c. 

We now show how to evaluate the GOE integral exactly. With the help of the identity ( D.l^ ), the 



GOE limit of j^ahp, as given in Eq. (5.8), can be written as 

-1 



I^goeP = 2 / dX I dXi I d\2 



1 + A 



M 



1 - A 1 



.Ai+As/ (Al + A2)2 7^ 
Further simplification results from observing that 

d fl \ _ 2(AiA2-A) 



(l + A)(AiA2-A) 



7^ 



dx\n 



7^2 



(C.14) 



(C.15) 



in which case 



L ^00 |.oo / 1 -|- A \ 1 — A 



+1 1 
dX- 

_i 1 + A 



1 - (A/ - 1) 



1 - A 
1 + A 



dXi / dA2 



1 


d 






xn 


^dX 




)] 


' 14 


M 


M+2 


1 


,Ai -f 


-A2y 




n 



(C.16) 



after an integration by parts. 

Next, we make the change of integration variables (Ai, A2) <—> (A, a), where 

Ai = l + i(A-l)(l + a) , 
A2 = l + i(A-l)(l-a) , 



the associated Jacobian being 



^(^ = i(A-l). 
5(A,a) 2V ) 



Then 



I^goeI^ 



+1 

da I dX 



1 - (M - 1 



1-A 



1 + A 



... /1 + A\*^+^ A-1 



dX 

1 Vl + A 



(C.17) 



(C.18) 



(C.19) 
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and we have 



7^ = (A - A)2 - 1(1 + A)(A - 1)2(1 - a^) 



(C.20) 



We are now in a position to implement the change of variables on the pair (A, A) given in Eq. ( |C.q ) (with 
Ai replaced by A there), followed by setting w — x{l ~ v). This leads to the representation 



I'S'goeP 



'^'^ \l + x)^f+^ 



where 



V 



1 .1 

Jo {I - v){x + 1 + v) + a^v^ ' 

da [ dv — — 5-^ . 

Jo {1 - v)ix + 1 + v) + a^v^ 



The a-integration can be performed, to yield the expressions 



eiix) 



dv ■ 



1 



: artan ■ 



v/(1 - v){x + l + v) y^{l-v){x + l + v) 

dv \ I : artan 



x + l + v ^(1 -v){x + l + v) 



(C.21) 



(C.22) 



(C.23) 



The integrals in Eqs. ( |C.23 ) are easily done for the special case of a; = 0, and give the results 
£i{0) = vr^/S, ^i(O) — ^2(0) = 1. For general values of x, we can proceed by making the change of inte- 
gration variable 



and introducing the quantity 



x/2 



f3^ 



Then we obtain 



eiix) 



du ■ 



p Vf~- u 
1 



1 + 2;/2 



: artan 



u- (i 



,2 ' 



The derivatives with respect to /3 have simple forms: 



1 — u u — P 

du , artan ■ 



d£i_ 

dh 
dp 



ln(l + /3) 
2/3 

1-/3 ' 
2/3 



ln(l + /3) 



where we have introduced 



(C.24) 



(C.25) 



(C.26) 



(C.27) 



e2{x)=£2{x)/{l + x/2) 



(C.28) 
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Now, after some integrations by parts, we can write 

^i(0)-^2(0) , 1 



1^, 



GOEl 



After noting that 



and 



M + 1 M + 1 (1 + xY^+^ 



dx ■ 



dx 



i'^ix) 



M- 1 



dx ■ 



I'^ix) 



M + 1 



(l + a;)M-i 2(M+l)yo {l + xy'+^ 







1 



2{M+1)Jq (3 \l + fi 



1-13- 



f3 

ln(l + f]) 



f3 



l-p y' 1 _ 1 d f l-(3 



M+l 



d_ 

'dp 



+ 



ln(l + (i) 



1-/3. 



ln(l + /3) 



(C.29) 



(C.30) 



(C.31) 



and using this to perform a further partial integration on the first integral on the RHS of Eq. (|C.29 ), it 
becomes clear that the two integrals in Eq. ( p. 29 ) cancel each other, leaving us with the desired result 



of Eq. (5.8) in the main text. 



Appendix D 

Preliminaries 

If we re-express TZ in terms of the variables {x,w,z) introduced in the previous appendix, and then make 
the further change of variables (w, z) ^ («', z'), where w = vx, z = 2xz' , we obtain the form 

n= .^f ,^ d:{z'), (D.i) 

(1 + vxY 

with 

Dl{z') = (z' + 1 + vxz'Y - 4(1 - v)z' . (D.2) 

Now, in the large-t asymptotic expansion, due to the exponential factor e^^*^ — e^®*^^ , the significant 
domain of the integration is confined to values xz' <ti 1, and hence vxz' <ti 1 (since < w < !)• Since also 
z' > 0, we can always make the replacement 1 — vxz' ~ 1 in our expression for D'^(z'), valid for large t. 
Thus, we have 

D^iz') iz' + lY-iil-v)z' = Dl{z'). (D.3) 

We introduce the notation Dy{z) = D^{z), and also note that one can write Dy{z) — {z — 1)^ + 4vz. 
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Likewise (in view of xz' <^ 1, vxz' ^1), we can make the following identifications and large-i re- 
placements: 

2x 



A1A2 — A — 



1 + vx 
2x 



[1 + (2 - v)xz' + z'' 
(1 + ^'), 



t»l 1+vx 
2x 

AA2 - Ai = [l + vxz'-z'] 

\ + vx 

t»i l + vx^'^ ^ ' 
2x 

A2-AA1 = 7:;— —[2v{l'^xz')^vx{l+vxz')-l + z'] 



(1 + vxY 
2x 



{2V-1 + VX + z) 
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t»i (1 + vxY 

We should note that, in these combinations of the A's, it is not justified to simply make the substitution 
A2 = 1 (i.e. z' — 0) in the asymptotic expansion, even when only leading order in t is sought, because it 
is not necessarily true that the significant integration region is confined to z' ^ 1. 



The following result will be used extensively in the subsequent discussion: Suppose that ]{x) is a 
function such that /(x 00) ^ foc{x) — Cx". Then, (i) if u < —1, we have 



/(te) 



(ii) if v > —1, we have 



dx 



f{tx) 



t" dx- 



foo{x) 



(1 + xY''+^ t»i 7o (1 + xY"""^ ' 



(iii) if = — 1, we have 



dx 



dxf{x) 



C 



C 



.7=1 ■' 



It is useful to note that with f{x 00) ~ Cx ^ and f{x) of the form 



fix) = / e-'-^g{z) 
Jo 



in which case we must have g{0) — 8C, Eq. (D.9) reads 
f{tx) 



,0 '^^ (1 + a;)^^+i t>>i 8s ""^ z Jo z 

where 7 denotes Euler's number. 



M ^ 

In8i + ln7-^- 
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It is also useful to define various functions for later convenience: 



1 , fl + z 



gi{z) = dv — 

Jo Dy{z) 4z 



= — In 



1- z 
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and we have gi(0) = 1, gi{z oo) ^ z 



and we have 52(0) — 1, 52(2 ^ 00) ^ z^^- FinaUy, 



4z 



93(2) 



dv 



2v 



1 



^?(^) 
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1 + z 



1 - z 



1 - z 
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Here we have 33(0) — 1, 53(2: — > 00) ^ z~^. Various integrals involving these functions can be established. 
These include 

/ dzgi{z)= [ dzg2{z)^'^, [ dz ^^-^ — - = 0. (D.15) 
Jo Jo ^ Jo ^ 

Also, 



dzgziz) = 1 



53(2) 



dz 

1 2 



(D.16) 



Large-t Expansion 



We isolate from the full expression (|5.2D, four distinct contributions to the 0{t ^) term in l^il 



U I \\ r d\^ /°°dA2e-2*('^-i) ( -ii^ 



A/+2 
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Ai = 



/+! i-oo foo / 1 
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.(l_A2 + A?-A^)e-2*(i-^^) 



1 + A 
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1 



(l + A)(AiA2-A) 



' d\ / dAi / dA2e-2*(^2-i) 

-1 ii ii VA1+A2; (1 + A)2 7^2 
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1 A1A2-A 



7^ 
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The quantity L\ contains the term of leading order in t ^ , i.e. half the GUE limit, and can be seen to be 
equivalent to the RHS of Eq. ( |C.l|) upon making use of the identity 



1 - A A2 + A2 - 2 



27^ 



1 



(l + A)(AiA2-A) 



7^ 
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1 + A (Ai+A2)2 (l + A)(Ai +A2)2 

We need to expand this term to first order is . The integrand of A\ comprises the terms whose 
i-dependence resides in the common factor of cxp[— 2/;(A2 — 1)]. This contribution is non-leading, and 
its large-i limit will yield a term of order 0(t^^). The integrand of B\ comprises the terms whose t- 
dependence resides in the common factor of 4t exp[— 2i(A2 — 1)] exp[— 2t(l — A^)]. This contribution is 
also non- leading, and its large-t limit will yield a term of order 0{ir^\ The term C\ would appear to be 
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of higher order in t ^. However, because of the singular nature of the factor TZ ^, naive power counting 
does not apply, and Ci does in fact contribute at order 0{t~^). There are four analogous contributions 



to the 0{t ^) term in the expression (5.4) for |S'2p, given by 
A2 



= U 



fl poo i>oo / -1 

d\j dX,J dA2e-2*(^^i) 



1 + A 



M+l 



A2 Ai — AA2 



+ A2; (1 + A)2 

- M+l ^ ^ 
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4t I dXJ dXi J 
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1 + A 

A/+1 



1 1 

(1 + A)2^ 



[(A2-AAi)(A2+A2-1)] e-2*(i-^^), 



dX / dAi / dA2e-2*(^2-^) . ^ ^ 

-1 Jl Jl V'^l + ^2 



1 + A \*^+^ 1-AAi(AiA2-A) 



1 + A 



7^2 
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One deals with Li and L2 most easily by combining them into L = _Li + £2- Then 

dX dxJ dA2e^2.(AM)[^ 1 + A \ A, 

-1 Jl Jl V'^i + A2 

/+1 roo POO / 1 I 

dX dXi dA2e-2*(^'-iM 
-1 Jl Jl V-^l + ^2 

\2iO('\ t \\l \ \\ \ f \ 1\2 



(1 + A)3 

1 + A XI 1 



(l + A)3 7^ 



(Ai - A)^ + 2(A2 + A)(Ai - A) - A(A2 - 1)^ 
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Next, we write L — I^gueP + which isolates the part that is zeroth order in < ^ . To extract (5i, 



which denotes the contribution of order 0(t ^), we must expand the integrand, save for e 



-4tz 



to 



first or second order in z = A2 — 1. We shall split the result of the expansion into three components, 
(5i — 6L' + 6L" + 6L'" , where SL' represents contribution from the first term on the RHS of Eq. ( D.20D , 



SL" is the contribution coming from the expression in square brackets in the second term on the RHS of 



Eq. ( D.20| ), and 6L'" denotes the rest. Then we have 

M+2 



6L' 
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dx ■ 
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At dx- 
Jo 



dXi 
2 



x 



dz e 
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t»l ./n (l+x)Af+2yo "'^'^'Jo D,{z) 



dz e 
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The contribution involving z'^ is non-leading, being of order 0{liit/t^), while in the rest we can set 
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Dy{z) ^ Dy{0) = 1 by virtue of Eq. to obtain 



SL" 



— I dx 



1 



t>i 8i Jo (1 + x) 



M+2 



dv (1 + vx) 
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16i 
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1 
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_M M + 1 

It is convenient to split 6L'" into a sum of two terms 5L'" = 5L'(' + SLi^' . For the first term, we write 



At I dX dXi / dz e""*^ , ^ ^ 
-1 Ji Jo VI + ^^1 



(l + A)3 7^ 



-2tz^ + 2z- (M + 2)- 



z 



[(Ai-A)2 + 2(l + A)(Ai-A)] 
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1 + Ai 

Due to the presence of the factor Ai — A in the numerator, it is safe here to set TZ ~z^o (Ai — A)^ for 



large t. This results in 



6L'I' 



, dX / dAi , 

4t./_i Ji Vl + Ai 



1 + A\'^+^ 1 



(1 + A)3 



2(1 + A) 
Ai - A 



M + 2' 



1 + Ai 



dx ■ 
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16t./n (l + x)^^+3 



dv [2(1 + x)- {M + 2)(1 + vx)] 



1 

16t 



-2 



1 



1 



M M+l_ 
The remaining contribution to SL'" is given by 
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dX / rfAi / dze 
1 1 



-4*wl + A\^^+^ 1 



1 + Ai 



7^ (Ai-A)2 



Using the relation 



(Ai-A)^ 



D^{z') - 1 + 2vxz' 
t»i Dy{z') 
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we find 



(5L!; 



2wa;z 
D,{z) 



(D.28) 



where the two components pick out each of the two terms in the square brackets, respectively. In <5i2b' 



we are allowed by Eq. (D^) to set Dy{z) ~ £'i,(0) — 1, and so obtain 



SL: 



2b 



dx ■ 



x + 2 



t»i 8t Jq (1 + x) 
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dv V 
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m 



1 



1 



M M + 1 



The component 5i2a can be expressed as 



SL'' 



dx ■ 



f{tx) 



{l + x)M+i 7o {1 + X)M+^ ' 
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where 



fix) = Ax dze-«^^[5i(z)-l] 



1 1 



In this case, Eq. (D.7) imphes that 



SL 



2a 



t>l t 



dx fix) — / dz 

■'^ ' 8t ' 



3 (8a;)2 
9i{z) - 1 



(D.31) 



(D.32) 



We see from Eq. (D.15) that (5L2a vanishes up to order 0{t ^), and consequently only SL^'^, survives in 
6L'^'. Collecting together Eqs. p.21\ ), ( p.23| ) ( p.25| ) and ( p.29| ), we obtain 



SL 

t>i 16t 



-2+A 



1 



M M + 1 
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To extract the large-t limit of Ai, we proceed as above by writing X2 = 1 + z and expanding the 
integrand to lowest order in z, except for e"'**^ and the combinations in Eqs. (D.4)-(D.5) which are kept 
intact. But first we re-express Ai as 

'■+^ roc roc ^ / 1 + A \*^+' 1 r(l + A)(AiA2-A) 



dX dxJ rfA2e-2*(^'-i) ( ^ 
"2A2(A2-AAi) 



+ A2 



(1 + A)3 
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1 



n 



- 1 



(D.34) 



It is now convenient to write Ai = ^J"'' + A']^^ + A^-^\ where A^"-* denotes the contribution to Ai, as 



given above, involving the factor 1/7?.". So, we have 
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dA 








/•OC 
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1 
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8t 





-itz 



1 + A 
1 + Ai 



Af+2 



(1 + AF 



Next, if we note that by Eqs. (D.4) and ( p.(^ ) we can make the replacement 

Ax 



2A2(A2-AAi)- (l + A)(AiA2-A) 



t»i (1 + vx)' 



[2{v- l) + vx\ 



then we have 
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+ 1 /"OO /'OO 
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1 Ji Jo 



-4.wl + A\^^^^ 1 1 
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(D.37) 
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The observation that 



2(A2 ~ AAi) 

7^ 



(D.38) 



(2) 

allows us to express A\ in a, form that is amenable to integration by parts. It is convenient to do so in 
order to avoid the emargence of singular integrals (convergent only as principal values) in the process of 
reducing to the large-i limit. Accordingly, we write 



• 4t 



(2a) 



d\ / dXi 
1 Ji Jo 

dX / dAi 

Ai - A 



M+2 
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r(AiA2-A) 



d ( 1 
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7^ 
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(2b) 



Let us first deal with A'^^^^ . We can write this component as 
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with 



/+1 /'Oc /• 
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-1 Ji Jo 
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Ai 1 

(l + A)2 7^ 



[(2 - v)x + 1] 



1 
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Now, the leading contribution to the term involving (2 — v), which is of order 0(t^^), is annihilated by 



the differential operator in Eq. ( D.40 ). In the rest, only the logarithmic term, whose presence follows 



from Eq. (D.9), survives. Thus, we obtain 

^(2b) 



1 



If d \ hit 
t»i 16 V dt J t 
_ 1 
16i ■ 



dze-**^"(?i(z) 



Finally, we have 



+ 1 



~U d\ dXi 



dz ( 



-4., 1 + A\-+^ A,-A 



,1 + Ai; (i + A)2 

With the help or Eq. (|D.27|) , this can be reduced to the form 



n (Ai - A)2 



^) ^ - 4t / dx 
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One should note that this expression for Al has a structure very similar to that of dLi^' in Eq. ( D.28 ) 
and an almost identical analysis leads to the result 

1 1 



t»i 16t M + 1 ■ 



(D.45) 
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Adding together all contributions to A^^"^ and combining thus with the results ( D.35| ) for A^"-* and ( D.37 ) 
for A^-^\ we see that 



A, ^ ± 1-1- 
t»i 16t V M 



(D.46) 



To extract the large-t limit of A2, we also split it into a sum of regular and singular parts, 



A2 = A'^2 ^ + where A^"^ denotes the contribution to A2 in Eq. (|D.19|) that involves 1/TZ". Then 



/+1 roc poc / 1 



1 + A 
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1 1 
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where the three components pick out each of the three terms in the braces, respectively. The first 
component yields 

r+l ,00 ,00 / 1 , \ \ M+l ^ _ ^2) _ _ ^) 
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having set Dy{z) ~ -Du(O) = 1 by virtue of Eq. (D.8). With the aid of Eq. (D.5), the second component 
gives us 



l(lb) 



-1 



-4tz 



d\ / dAi dz e 
-1 Ji Jo Vl + Ai 

00 

dx 
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The term involving z is non-leading, being of order 0(ln</t^), while in the rest we can set 
Dy{z) Dy{0) = 1 to obtain 



A 
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1 1 

t»i ~8t M 
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The contribution from the third component ^3^'^'' is of higher order 0{hit/t^), and hence discarded. On 



combining the three components, we obtain 
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Taking account of Eq. ( D.4 ), we have for A 
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We express this as 



where 



A. 



dx 
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(l + a;)^+i 
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Then, Eq. (D.7) apphed to f2{x) and Eq. (D.ll) appUed to fi{x), together with the integral identities 
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(D.16), imply that 
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The full result for A2 is obtained by adding Eqs. (D.51) and (D.55) 



The sum B = Bi + B2 can be simplified to yield 
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{l + x)M+i ' 

where we have used the fact that only values w = vx <^ 1 are significant for large t, and we have found 
it convenient to introduce 
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Application of Eq. ( |D.7[ ) to /i(a;) and /2(a;), followed by the x- integration, results in 
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These two contributions cancel out, and leave us with 
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according to Eq. ( D.ll ), having noted that in this application g{z) = 0(1 — z), with Q denoting the 
Heaviside step function. We note that B cancels the non-analytic term in ^2- 
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We combine C = Ci + C2 and note that 1 - + A? - A| = -Tl + 2Ai(Ai - AA2). Then we decom- 
pose C = C'-^-' + C'-^-' according to the power of TZ^^ that the contributions involve. Thus, we have 
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The term C*-^-* is clearly of higher order in l/t. To evaluate C^^\ we first integrate by parts in order to 

(2) 

avoid the appearance of singular integrals, as in A\ '. This yields 
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Making use of Eq. (D.27) gives us 
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Then, using Eq. (D.7) and performing the x-integration, we have 
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For the component C'^^'^-', we can write 
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The second term in the square brackets is non-leading, leaving us with 
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Hence, after appealing to Eq. ( D.7 ) and carrying out the a;-integration, we obtain 
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We see from Eqs. (D.65) and (D.69) that C^^-* vanishes up to order 0{t and therefore C provides no 
contribution to the final result. All the separate components can now be combined to give 

2 



SL + Y{A,+Bj+Cj) ^• 



M 



(D.70) 



which leads directly to Eqs. ( 5.10 ) and ( 5.11 ) in the main text 



Additional Remarks 



The double integrals which appear in Eqs. ( D.58 ), ( D.65 ) and ( D.69| ) can be easily performed by making 
the change of integration variables (u, z) h- s- (u, z) defined by 



Then we have 
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The results used in the text follow immediately. 



1 pOD pi pOO 

dv dz = I dv I dz , 
Jo Jo Jo {l + vz)3 

Dii~z) 



(1 + W5)2 ' 



(D.71) 



(D.72) 



Appendix E 



A numerical treatment of the integrals in Eqs. (5.2 5.4) cannot be implemented directly as their integrands 
exhibit an algebraic singularity on the surface 



TZ = X^ + Xl + Xl~ 2AA1A2 -1 = 



(E.l) 



which, of course, cancels out between numerator and denominator in the end, once all relevant contri- 
butions are taken into account. This apparent singularity can be eliminated altogether via the following 
change of integration variables: 



A = l~2q, 
Ai = 1 + 2qqi 
A2 = 1 + 29^2 



(E.2) 



Then q, qi, q2 span the ranges < qi,q2 < 00 and < q < 1- 
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Let us now consider the case t = and write 

\W= [ dqlM{q) . (E.3) 
Jo 

Then one can show, for example taking M = 1, that Ii{q) ^g^oo const, x so that there is an 

integrable singularity at g = 0. The divergence of luiq — > 0) does not disappear for larger values of M. 
It can be removed by the transformation g = l/(p + 1). The resulting integral is regular in p but displays 

a consequent slow fall-off as p ^ oo, going as p^^l"^ for M = 1. So, to achieve an accuracy of 10""*, the 
numerical cut-off in the p-integral would have to be taken as 10*. 



Nonetheless, we thus arrive at the expressions 



= Hdp r dq^ r dq2 (t—^ 

Jo Jo Jo \l+p + ( 



■qi+q2 



M-l 



exp[-8tq2{l+p + q2)/{p+l)^] 1 



+ +52)2 



1 + 91 + 92 + 



p 



29192 



p+ll+p + qi+q2 



p+1 



:(1 + £)p - (1 - S){qi - q2){l+P + 91 + 92)] 



(E.4) 



and 



JS^ = / dp dq, 

Jo Jo Jo 



dq2 



P 



1 + P + 91 + 92 



M-l 



exp[-8t(?2(l+P + 92)/(p+l)'] 1 



(l+p + qi+q2)^ 



(91 - 92) 



2p 



3 + ^ ((p + 1)(1 + 91 + 92) + 2^192) 



■{- 

+ (p ^"1)2 iiP + - 91 + 92) + 292) - (1 + p + 91 + 92)(1 +P + 92)' 
(1 - £)p 



p + 1 



+ 



p + 1 



1 + (b + + 91 + 92) + 29192) 



(l+p + 91 + 92) 



+ (p+l)2 [2(p+l)(l-9i+ 92)+ 49i + (9i -92)(l+P + 292)'] | • 



We have set 7?.p = (p + 1)7?. so that 



TZp = {p+ 1)(1 + qi+ 92)^ - 4p9i92 



(E.5) 



(E.6) 



and here £ = exp[— 8^ {p + 1)^]. In this form, the integration is completely free of singularities; but, for 
small t, the integrand decays slowly (with a power law) in all directions. Consequently, all three upper 
integration limits must be taken very large, and this makes accurate computation slow and difficult. 

To circumvent this, we introduce two compact coordinates by transforming to two angle variables a, 
w and one radial parameter r according to the definitions 

91 = \r{l - a){l+w) , 

92 = ir(l - a)(l - «;) , 

p = ra . (E.7) 
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The ranges of the new variables are given byO<r<oo, 0<a<l, —l<w<l. However, because of 
the very asymmetric form that the exponential now acquires, there exist both rapidly and very slowly 
decaying directions for the integrand. Indeed, the upper limit of the r-integral must still be taken 
extremely large in order to achieve sufficient accuracy for the slowly decaying parts. However, this means 
that the integration region also extends far along the rapidly decaying directions where the integrand 
is virtually zero; and, in fact, the integrand turns out to have support only in a tiny fraction of the 
entire integration volume. The integration routine is then unable to properly determine when the desired 
accuracy has been achieved. This is because, at some stage in the iteration, further subdivision of the 
grid is likely to create new points essentially only in the pervading insignificant region, and hence not 
alter the previous value of the integral to within the requisite accuracy. 



The easisest way to remedy this problem is by compactifying the x-integration through the variable 
change 



y = 



l + r 



(E., 



so that < y < 1. The transformations of Eq. (E.7,E.8) result in a simple bounded integration over a 
hyper-rectangular region and an integrand that has a significant value over an appreciable fraction of this 
region. The change of integration measure and region that corresponds to the variable transformation 
{p, qi, (72) (y, OL, w) can be summarized by the equation 

P \ 1 



00 poo poo / __ \ Af— 1 

dp dqi dq2 i — 
Jo Jo \l+P + qi+q2 



{l+p + qi+q2y 



1 



dy 



y 



M+l 



a-yy 



1 /•+! 

/ dw 

1 



da (1 — a)a^ 



(E.9) 



Then we can write for |S'jp (j — 1, 2) 



where 



Fi[y,a,w) 



and 

F2{y,a,w) 



1 



f dyv^'+^ f da{l-a)a^"-^ T 
Jo Jo J -I 



■ exp I —2t 



y{l - - w) 
l-y{l-a) 



7^„ 



y{l - a){l - w) 



1-2/(1 -a) 



1-y 



ly^a(l~af(l-w'^) 
1 - - o 



2 l-y(l-a) 



2/(1 - y(l - ")) 



Fj{y,a,w) 



[{1 + £){! ~ y)a ~ (1 ~ £){! ~ a)w] 



+ 4t 



y(l ~ y)a , / y[\ — a){\ — w) 



\\-y{\-a)Y 



1 - y(l - a) 



y{l — a)w 

i-y(i-a) 



2ya 



(E.IO) 



(E.ll) 



3 + ^ ((1 - y) + y'a{l -a) + \y\l af{l - w^)) 
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2ya [l — y — y^a{l — a)w) — (1 — 2/(1 — a) 



2(l-y)(l-y-2/2a(l-a)«;) 



At 1 

+ 



with 



1 - 2/ [1 - y(l - a)]2 
+ 2/(1 — a)w (1 — 2/(1 — a)tt;)^ 



Tly = 2/'(l-2/)(l-a)' + (l-2/)[l-2/(l-a)][l-2/ + 22/(l-a)]+2/Ml-a)^"' , 
2/(1 - y)a 



(E.12) 



f = exp < —8t 



[l-2/(l-a)]2 



_ 2/(1 

In the numerical analysis, we represent as the integral 

jr(4te) = 2 [ due"^*^" (E.14) 

^0 



when the argument is small. It is also useful to note the relation 

1 {i-yf 



(E.15) 



In the parameterization above, a mild algebraic singularity appears along the line y = 1, w = 0. This can 
be easily dealt with by setting the upper limit of the 2/-integration to 1 — 6 with 6 an arbitrarily small 
positive number. 
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Figure captions 



Figure 1: Representation of a typical lead-stadium configuration. The size of the stadium L is large 
compared with the clastic mean free path ensuring that the electron motion is ballistic. A total of M 
external channels couple to the stadium through the two leads. 

Figure 2: Schematic depiction of the zero- field dip in the average magnetoconductance associated with 
weak localization. The curve interpolates between values associated with pure GOE and GUE for low 
and high magnetic fields B, respectively. 

Figure 3a: Graph of the A3 statistic versus energy interval L for flux angle ^/t/io = 0. The solid line 



corresponds to setting t = in Eq. (3.4) (i.e. the GOE result). The dashed line is the GUE result. 



Figure 3b: Graph of the A3 statistic versus energy interval L for flux angle (/)/(/'o = 0.1. The solid line 



shows the best fit to the data points furnished by Eq. (3.4), obtained by adjusting t — 0.15 The upper 



and lower dashed lines correspond to the GOE and GUE results, respectively. 

Figure 3c: Graph of the A3 statistic versus energy interval L for flux angle 4'/4>o = 0.2. The solid line 



shows the best fit to the data points furnished by Eq. (3.4), obtained by adjusting t = 0.6 The upper and 



lower dashed lines correspond to the GOE and GUE results, respectively. 



Figure 4: The correlator K{t) (defined in Eq. ( |3.9| )) as a function of the time t computed as a numerical 
average over phase space, taking 10^ trajectories. The t ^ 00 asymptote is k ~ 0.06. 

Figure 5: Graph of the reciprocal of the weak-localization term (5(/^^ versus the magnetic-field parameter 
t for M = 1, 2, 4, 6, 8, 10 on a large scale (0 < t < 20). 

Figure 6: Graph of the reciprocal of the weak-localization term 5g~^ versus the magnetic-field parameter 
t for M = 1, 2, 4, 6, 8, 10 on a smaller scale (0 < t < 4). 

Figure 7: Plot of the weak-localization term 5g as a, function of ^/t for the case A/ = 1, compared with 
the corresponding Lorentzians obtained from fitting to the large--\/i tail of the M ~ 1 curve, and to its 
FWHM. 



48 



L 




Figure 1 




Figure 2 



o 

CD 



J O 



J O 



J o 

CO 



J o 

(>2 



_ O 



-I 1 1 1— 



-I 1 1 u 



CD 
O 



o 



o 



o 
o 



Reciprocal of g(GUE,M) - g(t,M) vs t 




0.5 1 1.5 2 2.5 3 3.5 4 



